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Spectral singularities are spectral points that spoil the completeness of the eigenfunctions of cer- 
tain non-Hermitian Hamiltonian operators. We identify spectral singularities of complex scattering 
potentials with the real energies at which the reflection and transmission coefficients tend to infinity, 
i.e., they correspond to resonances having a zero width. We show that a wave guide modeled using 
such a potential operates like a resonator at the frequencies of spectral singularities. As a concrete 
example, we explore the spectral singularities of an imaginary PT-symmetric barrier potential and 
demonstrate the above resonance phenomenon for a certain electromagnetic waveguide. 
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I. INTRODUCTION 



II. SPECTRAL SINGULARITIES 



Complex PT-symmetric potentials [l| that have a real 
spectrum are interesting, because one can restore the 
Hermiticity of the corresponding Hamiltonian and up- 
hold unitarity by modifying the inner product of the 
Hilbcrt space [2hJ]. It is usually believed that one can 
similarly treat every non-Hermitian Hamiltonian H that 
has a real and discrete spectrum. This is actually true 
provided that H has a complete set of eigenvectors 0- 
For the cases that the spectrum is discrete the lack of 
completeness is associated with the presence of excep- 
tional points. These correspond to situations where two 
or more eigenvalues together with their eigenvectors co- 
alesce. This phenomenon is known to have physically 
observable consequences It also plays an important 
role in the study of open quantum systems particularly 
in relation with the resonance states @, 0] • For the cases 
that the spectrum has a continuous part, there is an- 
other mathematical obstruction for the completeness of 
the eigenvectors called a "spectral singularity" [22[ • The 
purpose of the present article is to describe the physical 
meaning and a possible practical application of spectral 
singularities. 

Spectral singularities were discovered by Naimark Q 
and subsequently studied by mathematicians in the 
1950's and 1960's @. The mechanism by which spectral 
singularities spoil the completeness of the eigenfunctions 
and their difference with exceptional points are discussed 
in 0. 

Spectral singularities of complex 'PT-symmetric and 
non-PT-symmetric scattering potentials have been stud- 
ied in pjj HH • In this article we shall examine the spec- 
tral singularities of the imaginary potential [H, EH : 



il for - a < X < 0, 
-il for < X < a, 
otherwise, 



(1) 



with a € M. + and 3 e M — {0}, that has applications in 
modeling certain electromagnetic waveguides (l2j . 



Consider a complex scattering potential v(x) that de- 
cays rapidly as — > 00 [23|. Suppose that the con- 
tinuous spectrum of the Hamiltonian H = — + v(x) 
is [0,oo), and for each k e M+ let tpk±( x ) denote the 
solutions of the eigenvalue equation Hip(x) = k 2 ip(x) 
satisfying the asymptotic boundary conditions: 



ipk±(x) -> e ±tkx as 14 ±00, 



(2) 



i.e., the Jost solutions. A spectral singularity of H (or 
v) is a point k 2 of the continuous spectrum of H such 
that the ipk+± are linearly-dependent, i.e., they have a 
vanishing Wronskian, V'fc.+V'fc.- ~ V'fc*-^ + = 0, (Toj |- 

Clearly the continuous spectrum of H is doubly- 
degenerate. To make this explicit, we use ip 3 . with k £ IR + 
and q £ {1, 2} to denote a general solution of the eigen- 
value equation Hip(x) = k 2 ip(x). Furthermore, because 
v(x) — > as x — > ±00, we have 



lkx + si. 



as 1-) ±00, 



(3) 



where A± and B± are complex coefficients. A quantity 

of interest is the transfer matrix M(fc) that is defined by 

A B \ ( A B \ 

J ) = M(fc) ( ^ ) . Among its useful properties 



are the identity det M(fc) — 1 and the following theorem. 

Theorem 1: k~. € M + is a spectral singularity of H 
if and only if either —k+ or fc* is a real zero of the 
entry M 22 (fe) o/M(fc), 0. 



Next, consider the left- and right-going scattering so- 



lutions of Hi\){x) = k 2 ip(x) that we denote by ij) l k and 
respectively. They satisfy 



N r 



Jkx + / .', 

NiT l e ikx 

NrT r p -ikx 

+ R r e 



ikx\ 



-ikx 



r a ikx\ 



as x 
as x 



as x 



—00, 
+00, 

> —00, 

> +00, 



(4) 
(5) 
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where Ni, N r , R l , R r ,T l and T r are complex coefficients. 
N[,N r are normalization constants, |i?'| 2 , \R r \ 2 are the 
left and right reflection coefficients, and |T / | 2 ,|T r | 2 are 
the left and right transmission coefficients, respectively. 
Comparing ((4]) and ([5]) with ([2]), we see that tp l k and 
ip k are respectively proportional to the Jost solutions 
ipk+ and ij)k- ■ Therefore, at a spectral singularity, 
k 2 , the scattering solutions ip l k and ip^ become linearly- 
dependent. In view of (HJ and ([5]), this is possible only 
if R l ,R r ,T l and T r tend to infinity as k -> fc*. The 
converse of this statement is also true: 

Theorem 2: fc 2 € R + is a spectral singularity of H if 
and only if the left and right reflection and transmis- 
sion coefficients tend to infinity as k — > fc* or k — > — fc*. 

The following is an explicit proof of this theorem. 

Comparing (J4j) and §5§ with J3|) , we can determine the 
coefficients A± and £>^_ for ijj k and and use them to 
express R l , R r , T z and T r in terms of the entries of the 
transfer matrix M(fc). This yields 



T l = 1/M 22 (fc), 
T r = 1/M 22 (k), 



R l 
R r 



-M 21 (k)/M 22 (k), (6) 
M 12 {k)/M 22 (k), (7) 



where we have employed det M(fc) = 1. As seen from © 
and ([7]), at a spectral singularity, where M 22 (k) vanishes, 
R l ,R r ,T l and T r diverge. The converse holds because 
Mi 2 and M21 are entire functions (lacking singularities). 

Another curious consequences of ffib and (J7J , is the 
identity: T l — T r . This is derived in [14| using a different 
approach, but is usually overlooked. See for example [l5| . 

Next, we examine the S'-matrix of the system: S = 

T i ), M- In vi cw of ©, ©, and dctM(fc) = 



1, the eigenvalues of S 



are given by s± = (1 ± 
M\i{k)M 22 {k))/M 22 (k). At a spectral singular- 
ity s + diverges while s_ Mn(k)/2. This suggests 
identifying spectral singularities with certain type of res- 
onances. Indeed, in view of Theorem 2 and Siegert's 
definition of resonance states [lj|, they correspond to 
resonances with a vanishing width (real energy) . 



III. VT- SYMMETRIC BARRIER POTENTIAL 



Consider the Hamiltonian operator H = — j^+v at5 (x) 
with v a ^(x) given by (P). Because v a ^(x) = for |x| > a, 
the results of Section I I apply to v a>} . The determination 



of the eigenfunctions [13|, [ljj of H and the correspond- 
ing transfer matrix M(fc) is a straightforward calculation. 
Here we report the result of the calculation of M 22 (fc): 



M 22 (k) = e 2iak [h(k) - if 2 (k)]/Vl+7, (8) 
where f\ and f 2 are real- valued functions given by 
fx(k) = y/1 + y 2 |cos(afcw)| 2 - |sin(ofcw)| 2 , (9) 



f 2 (k) = di y/l + iy(2 - iy) sin(akw) cos(akw*) , (10) 



y --=t/k 2 



w 



— iy, and 5ft means "real part of" 



According to Theorem 1 and Eq. ((U, k 2 E M. + is a 
spectral singularity of v a . } if and only if fi(k) = and 
f 2 (k) = 0. If we insert © and (flU|) in these equations 
and divide their both sides by | cos(afciu)| 2 , we find 



|tan(ofciy) 



tan(akw) = — 



VT=iy{2 + iy) 



VTTiy(2 - iy) 



(11) 

tan(afcw)*. (12) 



Now, we multiply both sides of (fT2")) by tan(afciu), use 
(HTJ, cos(26») = (1- tan 2 6> )/(l+tan 2 9) and w = y/l -iy, 
to obtain cos(2afci/l - iy) = -(1 + 4y~ 2 ) + 2iy~ 1 . This 
equation is equivalent to 

cos r cosh q 
sin r sinh q 

where 



-(1- 
22T 1 . 



4?T 2 ), 



ak 



ak 



1-1 sgn(y), 



y 2 + 1 + 1 



(13) 
(14) 



(15) 



(16) 



sgn(y) denotes the sign of y, and we have employed 
the identities sin(^2 
and cos(- tan ~ 



) = sgn(yWi[l-(t/ 2 + l)-l/2] 



2 [l + (ya + l)-i/a]. 
Next, we solve for y^ 1 in (fM)) . substitute the result- 



ing expression in (|13p . and use the identities sinh g = 
cosh 2 q — 1 and cos 2 r = 1 — sin 2 r to obtain a quadratic 
equation for coshg with solutions 

1 



coshg = ^(— 1 ± ^/2cos(2r) — l)cotr cscr. (17) 

To ensure that the right-hand side of this equation is 
real, we must have cos(2r) > i. Furthermore according 
to (fT3"l). cos(r) < 0. These imply 

7T 

\r — (2n + l)n\ < — , for some integer n. (18) 

Under this condition the right-hand side of (|17l) is 
greater than 1. Hence, q = ±q±(r), where q±(r) := 

-1 



cosh 



If we set 



cotr cscr(± v /2cos(2r) - 1 - l)/2 

q = ±q±(r) in ()14j) and solve for y, we find y = 
± sgn(sinr)y±(r), where y± := 2\ sin r sinh q±(r)\ . 
Inserting this expression for y in (jl5p and (|16p and 
solving for q give q = ±q±(r), where q±(r) := 



rsgn(sinr) 



y / V±(r) 2 + l-l 



The spectral singularities cor- 



y/y±(r) 2 + l + l' 

respond to the values of r for which q± (r) = g± (r) . These 
are transcendental equations admitting simple numerical 
treatments. It turns out that q+(r) = q+(r) does not 
have a real solution fulfilling ([T8]). while <7-(r) = <?-(r) 
has two solutions ±r„ for each choice of n in (|18p . Ta- 
ble 1 lists the numerical values of r n for various choices of 
n. It turns out that r n > and r_„ = — r n +i for n > 0. 
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n 


r n 


y n 


ak n 


a 2 in 





2.64390700 


1.82765566 


1.06468255 


2.07173713 


1 


9.11655393 


0.71364271 


4.31823693 


13.3074170 


2 


15.4804556 


0.49008727 


7.52928304 


27.7830976 


10 


65.8884385 


0.17167639 


32.8243878 


184.971084 


100 


631.445619 


0.02901727 


315.689592 


2891.85852 



TABLE I: r„, y n , k n , and 3„ are respectively the numerical 
values of r,y~,k and 3 that correspond to spectral singular- 
ities. For n > 0, k n and j n are increasing functions of n. 
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FIG. 1: Cross section of a waveguide with gain (+) and loss 
(-) regions in the x-z plane. Arrows labeled by /, R, and T 
represent the incident, reflected, and transmitted waves. 



Next, we insert y = ± sgn(sin r) y_ (r) in (|16p and 
use the identity a 2 3 = {ak) 2 y to obtain k = g(r) 
and a 2 j = ztg(r) 2 sgn(sinr) y_(r), where <?(r) := 

r / 2(^y^(r) 2 +T + 1). Setting r = ±r n in these re- 
lations gives the values (ak n and a 2 3„) of a/c and o 2 3 that 
are associated with spectral singularities. We list some 
of these values in Table Q] Because ak and a 2 j are odd 
functions of r and r_„ = —r n +i f° r n > 0, we have 
k- n = —k n +i and 3_„ = — 3 n +i- According to Tabic HI 
the smallest values of a|fc| and o 2 |3| for which a spec- 
tral singularity occurs are respectively afcrj ~ 1-06 and 
a 2 3o ~ 2.07. Using more accurate values for ak n and 
a 2 3 ra that we do not report here, we have checked that 
\M 22 (k n )\ < 1CT 9 for |n| < 20 and \n\ = 10 2 , 10 3 , 10 4 . 

For the system we considered in this section, each value 
of a 2 3 can support at most one spectral singularity (either 
the latter does not exist or it exists for a single energy 
value) . 

IV. A VT- SYMMETRIC WAVEGUIDE 

Consider a rectangular waveguide with perfectly con- 
ducting walls that is aligned along the z-axis and has 
height 2/3 as depicted in Figure [TJ Suppose that the re- 
gion \z\ < a inside the waveguide is filled with an atomic 
gas, and a laser beam shining along the y-direction in 
the region — a < z < is used to excite the resonant 
atoms and produce a population inversion. In this way 
— a < z < and < z < a serve as gain and loss re- 
gions respectively, and the relative permittivity at the 

resonance frequency takes the form e(z) = 1 + — E ^p—— - 
for \z\ < a and e(z) = 1 for \z\ > a, where io,io p and 
(5 are respectively the frequency of the wave, plasma fre- 
quency, and the damping constant, Alternatively, 
e(z) = 1 — v a ^/ UJ (z) where s := to 2 / (25). In [lU, the au- 
thors use an approximation scheme to reduce Maxwell's 
equations for this system to the Schrodingcr equation for 
the barrier potential (p}. Here we offer an exact treat- 
ment to examine singularities of the reflection and trans- 
mission coefficients for this waveguide. 

Let i,j,k be the unit vectors along the x-, y- and 
z-axes, & := us/c, m € Z + , &, n := nm/(2/3), 
Xm{x)) := sinl&mix + /?)], and k := ^/ & 2 - & 2 n . Then 



E(r,t) = $l(e- iut [-i(JXm(x)4>(z)]j) , and B(r,t) = 

(e~ zut [xm{x)(f>' (z)i — i arc transverse 

electric (TE) waves satisfying the boundary conditions 
for the waveguide and solving Maxwell's equations pro- 
vided that 

<j>"{z) + [^E{z)-^ n ]4>{z)=Q, (19) 

and <f> and <p' are continuous functions on the z-axis 
[Hj]. For \&\ > the solution of has the form 
([3| with z and k playing the roles of x and k respec- 
tively. This allows us to define a transfer matrix M(k) 
for this system and introduce the right and left transmis- 
sion and reflection amplitudes, T l r and R Lr , associated 
with (HU). These satisfy © and ©, and diverge when- 
ever M 22 (k) = 0. 

It is not difficult to see that the right and left reflec- 
tion and transmission amplitudes for the propagating TE 
wave coincide with T l ' r and R l ' r , respectively. Therefore, 
if we can tune the frequency uj of the incoming wave to 
the frequency of a spectral singularity, then the am- 
plitude of the wave will diverge as uj — > cj*. In practice, 
this means that sending in a wave of frequency uj « 
will induce outgoing (transmitted and reflected) waves of 
considerably enhanced amplitude. The waveguide then 
uses a part of the energy of the laser beam to produce 
and emit a more intensive electromagnetic wave. Note 
that this effect is fundamentally different from the res- 
onance effects associated with exciting resonance modes 
of a cavity resonator. Unlike the latter that has a geo- 
metric origin, the spectral singularity-related resonance 
effect relics on the existence of a localized region with a 
complex permittivity (a complex scattering potential). 

The calculation of the transfer matrix M(re) defined 
by (|19[) is analogous to that of the "PT-symmetric bar- 
rier potential. In fact, M 22 (k) takes the form ([5]) pro- 
vided that we set: k = n, a = a, and y = s&/(ck 2 ). 
In particular, we can determine the values of uj and 
s for which the above resonance phenomenon occurs 
by setting k — k n and s&/(ck 2 ) = ±y n . This yields 
w = Lo n , m and s = s n ,m where w n>m := c^J k 2 + $? m and 
S n ,m ■= ±ck 2 t y n / '^Jk 2 +& 2 n = ±c 2 3„/w„, m . For m = 1, 
huip = 0.2 eV, hS = 1.25 eV, we attain the spectral sin- 
gularity with n = for fko = Hluq i = 5 eV, a ~ 1004 nm 
and (3 « 62 nm. Figure [5] shows the graphs of the log- 
arithm of the transmission and reflection coefficients as 
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FIG. 2: Graphs of log 10 (|T r '' ! | 2 ) (full blue curves), log 10 (|7?' | 2 ) 
(dotted red curves), and log 10 (|i? r | 2 ) (dashed green curves) 
as a function of oj/u>o,i, for m = 1, n = 0, hwo,i = 5 eV, 
tiwp = 0.2 eV, M = 1.25 eV. For the figure on the left (right) 
a = 1004.17 nm, j3 = 62.0464 nm (a = 1004 nm, j3 = 62 nm). 



functions of w/u>o,i for this case. The location and hight 
of the pick representing the spectral singularity arc highly 
sensitive to the values of a and /3. 

We close this section by noting that similar VT- 
symmetric waveguides have been considered in [IH, . 
These differ from the one we studied in that in our case 
the permitivity changes along the direction of the propa- 
gation of the wave (z-axis). This is essential for realizing 
the spectral singularity- related resonance effect. 

V. CONCLUDING REMARKS 

In this article, we offered for the first time a sim- 
ple physical interpretation for the spectral singularities 
of complex scattering potentials. In the framework of 
pseudo-Hermitian quantum mechanics [|[ , where one de- 



fines unitary quantum systems with a non-Hermitian 
Hamiltonian by modifying the inner product of the 
Hilbcrt space, the presence of spectral singularities is an 
unsurmountable obstacle [Ioj |. In contrast, in the stan- 
dard applications of non-Hermitian Hamiltonians, spec- 
tral singularities are interesting objects to study, because 
they correspond to scattering states (with real energy) 
that nevertheless behave like resonant states. 

We explored the spectral singularities of a VT- 
symmetric potential v a<2 that admits a realization in the 
form of a waveguide. We obtained the values of the phys- 
ical parameters of the waveguide and the propagating TE 
wave for which the system displays the resonance behav- 
ior associated with the spectral singularities of v a:} . 

Our results call for a more extensive investigation of 
the spectral singularities of complex scattering potentials 
that can be realized experimentally. This should provide 
means for the observation of the resonance effect that 
is foreseen in this article. Another line of research is to 
explore the spectral singularities of complex periodic po- 
tentials [2(| ■ A more basic problem is to study the conse- 
quences of spectral singularities for the implementation 
of quantum scattering theory. Similarly to exceptional 
points, presence of spectral singularities leads to sub- 
tleties associated with the existence of an appropriate res- 
olution of identity [2l[ . This problem may be avoided for 
scattering wave packets obtained by superposing eigen- 
functions of the Hamiltonian. A general treatment of this 
problem requires a separate investigation. 
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